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1 Introduction 

X. 

In this paper we consider a nonisothermal multi-phase transition process occurring in a 
bounded container Q C K , iV G N, with Lipschitz boundary dQ. The state variables 
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describing the evolution of the system are the absolute temperature 9 > and the vectorial 
order parameter \ £ Mr, d £ N. Following the idea that was already described in the pio- 
neering papers [26] and [I] , but which has been only recently analyzed in a more systematic 
way (cf., e.g., [I] [2], [5] [6] , [8] [21], [21]), we take into account long range interactions 
between particles. Then the model equations resulting from the energy and entropy balance 
relations have the form 

(1.1) (e(9, X ))t + (A(x) + 0<f(x)) t + b[x\Xt ~ div (k(6, X )V0) = in Q M := Q x (0, +oo), 

(1.2) n{9) X t + A'(x) + b[ X ] + ((3 + 0)dip{ X ) + Oa'ix) + e x (0, x) - 9s x (9, X ) 3 in Q,, 

(1.3) fc(0,x)V0-n + 7(0-0 r ) = on Soo := dVt x (0, oo), 

(1.4) 9(-,0) = 9 , x(-,0) = xo infi, 

where n denotes the outward normal vector to dQ, and (11.21) has to be understood as an 
inclusion in M d , where dip is a possibly multivalued subdifferential of a general proper, 
convex, and lower semicontinuous function </? : IR d — >■ M U {+00} . The physical meaning of 
the functions e, s, A, a and of the positive constant (3 is explained in (jl.l6p - (jl.l9p . while 
b[x] (whose explicit form will be given) represents the nonlocal operator acting on x- With 
an abuse of notation we have used the symbols A' , a' , e x , s x for gradient vectors in M d 
and omitted the scalar product symbol between M. d vectors (like b[x] and Xt i n fll-ip ) in 
order not to overburden the presentation. The function \i in (11.21) represents the (bounded 
away from 0) mobility of the system, while 7 denotes the heat transfer coefficient through 
the boundary dQ. The external temperature 9r is a sufficiently regular boundary datum 
on Sqo , and # , Xo are supposed to be two given initial configurations. 

The main novelty here is to consider a multi-phase nonlocal phase field system in 
the case when the specific heat cy(0,x) = dee(9,x) and the heat conductivity k(9,x) are 
not constant and depend on both the variables 9 and x- Suitable regularity and growth 
conditions will be specified in the following section. 

Let us only note that many typical expressions for cy in a two-phase system (i.e. in 
case d = 1) can be included in our analysis. In the solid-liquid system mentioned above, for 
example, we may have different values c v {9) in the solid and c v {9) in the liquid phase, hence, 
we may define cy(9, x) — c i/(#) + x( c y(^) — c v(@)) [23 Section IV.4]). The value of x can 
be kept between and 1 by setting ip = i] 01 ] (the indicator function of [0, 1] ). The physically 
meaningful case in which the behaviour of c v and c v are powers of 9 ( ~ 9 a , a > 1 ) near 
zero and bounded functions for large 9 's can be covered by our analysis. Regarding the heat 
conductivity k, typical expressions of the type k(9,x) = Ki(9)x + K 2 (9)(l — x) , in case of 
a two phase system with x [0, 1] , for quite general functions K\ and K2 , are also allowed 
here. 

The main goal of this paper is to study the global existence of solutions to system 
(jl.lj) — (jl.4p . coupling a suitable variational formulation of the semilinear parabolic partial 
differential equation (11.11) for 9 to the integrodifferential inclusion (11.21) for x- We also 
prove some uniform in time upper bound for the absolute temperature of the system (see 
Theorem 12.21 below). Uniqueness of solutions is obtained under additional assumptions, in 
particular, in case that the heat conductivity k in (II .ip depends only on 9 and not on x- 

Before entering into the mathematical discussion of the problem, let us give a brief 
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derivation of the system fll.lj) — (jl.4p . emphasizing, in particular, the differences between local 
and nonlocal models. 

We assume here that the multi-phase transition process can be completely described 
by the evolution of the state variables 9(x, t) > 0, which represents the absolute temperature 
of the system, and the order parameter x( x ,t), which here is a vector in M. d . We fix some 
constant reference temperature 9 C , which will be assumed to be equal to 1, for simplicity. 

Inspired by the nonlocal Cahn-Hilliard model studied by Gajewski in [9], we consider 
the following nonlocal specific free energy 

F[9, X )=fo(0,x) + B[x], 

where B is a potential that accounts for long range interaction between particles. More 
specifically, given a bounded, symmetric kernel k : Q x Q — y IR and an even smooth 
function G : M. d — > M, we choose 

(1.5) B\x](x,t) := [ K(x,y)G(x(x,t)- X (y,t))dy. 

Jn 

Note that the local potential (u/2)\Vx\ 2 used often in the literature, see [25] and the refer- 
ences therein, can be obtained as a formal limit as n — > oo from the nonlocal one with the 
choice G(r]) = \rj\ 2 /2, n(x,y) = n N+2 K(\n(x — y)\ 2 ) , where k is a nonnegative function with 
support in [0, 1] and v = 1/N J RN k(\z\ 2 )\z\ 2 dz . This follows from the formula 

/ n N+2 K(\n(x - y)\ 2 ) \ x {x) - x{y)? dy = [ £(M 2 ) 
Jn Jn„(x) 

^ / K{\z\ 2 )(y X {x),z) 2 dz = u\V X {x)\ 2 
Jrn 

for a sufficiently regular x, where we denote fi n (x) = n(Q — x) . We have used the iden- 
tity j RN k(\z\ 2 ) (v , z) 2 dz = 1/N J RN k(\z\ 2 )\z\ 2 dz for every unit vector v G M. N (cf. the 
Introduction of [TH] for further details on this topic). 

Let E and S be the total energy and entropy densities, respectively. The process 
is governed by the internal energy and entropy balance relations over an arbitrary control 
volume CcO, 

(1.6) -A / E ( 6 > X)dx+ [ (q, n) ds(x) = tf(ft') , 

at Jn' J on' 



x( x + £)~ x( x ) 



dz 



(1.7) A I S(6, X )dx+ ! (?n) ds(x) > 0, 

where q is the heat flux vector, n is the unit outward normal to dQ' , and \I/(f2') is the 
energy exchange through the boundary of Q' due to the nonlocal interactions. Since B[x] 
is a potential field, it does not contribute to the entropy production in the Clausius-Duhem 
inequality ( II. 7p . 
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The local form of the entropy balance reads 

^(0, X ) + divq-^^>O, 

and it has to be understood in the regularity context of Theorem [272] below. This is certainly 
satisfied if 

(q,V0> <0, 
9S t (9, X ) + divq>0. 

Assuming 9 > and a suitable regularity with respect to time (this will have to be justified 
in the next sections), we obtain from fll.6p that 



(1. 



Differentiating the identities F = E — 9S = fo + B[x] with respect to t, we obtain 

(1.9) F t = E t - 9S t - 9 t S = d e fo9 t + d x f oX t + B[ x ] t , 

where d x fo stands for an element of Clarke's partial sub differential of fo with respect to 
X G M. d , and dgf is the partial derivative of f with respect to 9 el. Consequently, 

(1.10) S = -def = s , E = e + B[ X ], f = e -9s , 
and inequality (I1.8P reads 

(l.H) / (d x f Xt + B\x]t)dx<*(tt). 

The nonlocal interaction takes place only inside the domain Q, hence = 0. A canonical 

way to satisfy these conditions independently of the evolution of x consists in choosing the 
order parameter dynamics in the form 

(1.12) f,{9) X te-D x T[9,x) 
with a factor /i(6 l ) > 0, where we denote 

T[9,x] = [ F[9,x] dx 
Jn 

and D X T stands for the Clarke subdifferential of J 7 with respect to the variable x G 
L 2 (n;M. d ). The inclusion sign in fll,12p accounts for the fact that fo(9, X ) includes terms 
that are possibly not Frechet differentiable. Condition fll,12p is based on the assumption that 
the system tends to move towards local minima of the free energy with a speed proportional 
to \/fi{9). Denoting 

(1-13) b[ X ]{x,t) := 2 / K(x,y)G'(x(x,t)-x(y,t)) dy , 
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where again, with an abuse of notation, G' stands for the (^-component vector VG, we 
see that the inequality (II. lip holds without prescribing any relationship between fi(8) and 
B[x], provided that we choose ^(SY) in (11.61) as 

(1.14) [ (-b[ X ]xt + B[x}t) dx. 

The differential form of the energy balance (II .6p then reads 

(1.15) Et + divq = -b[ X ]xt + B[x}t. 

The specific heat cy(8,x) is the only thermodynamic state function, which can be 
identified from the measurements, while the local internal energy and entropy densities are 
computed from the formulas 

(1.16) e {e,x) = e {0, X ) + e(6,x), e(8, X ) = [ c v (r, X )dr, 



1.17) so(9,x) = s (0,x) + s(9,x), s(0, X ) = I dr 



o ' 

where e (0, %), s (0, x) are i n fad "integration constants", which we choose as 

(1-18) eo(0,x) = A(x) + Mx), 

(1-19) S (0,x) = -<T(X)-<P{X), 

where ip : M. d — > K U {+00} is proper, convex, and lower semicontinuous, the functions A 
and a are sufficiently regular on T>(tp) , and the parameter f3 is a positive constant. 

Then, the free energy functional F has the form 

(1.20) F[9, X ] = f(9,x) + Hx) + B[ X } + W + 9Mx) + 9(T(x), 

where f(6, X ) = e(9, X ) - 9s(9, x) ■ 

Using (11.201) . we rewrite the phase dynamics (11.121) as 

(1.21) fM(6) X t + A'(x) + b[ X ] + W + 8)d V ( X ) + 8a' (x) + e x (8, X ) - 9s x (8, X ) 3 , 
while the internal energy balance (11.151) can be reformulated as 

(1.22) (e(8, X ))t + (A(x) + Mx)\ + b[x]Xt ~ div (k(8, X )V8) = 0. 

We now show that in the energy conserved case, that is, if we assume no-flux boundary 
conditions (7 = in (11.31) ). the phase transition model with a nonlocal interaction potential 
is compatible with the Ottinger-Grmela GENERIC formalism (abbreviation for "General 
Equation for the Non-Equilibrium Reversible-Irreversible Coupling"), see 
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Set 



1.23) £[6,x]{t) = [ E(0, X )(x,t)dx 

Jn 

1.24) S[e, X ](t) = [ S(6, X ){x,t)dx, 

Jn 

1.25) B[ X ](t) = f B\x](x,t)dx. 

Jn 



We show that there exists a symmetric positive semidefinite matrix M[#,x] such that 
(1.26) M[9, X ] 



DgE[9, X ]\ = (0 

D x s[e, x ]J V° 

and such that the system fll.2ip - fll.22p has the form 



It suffices to choose (we omit the arguments of the state functions for simplicity) 

0/ lmi2 m 22 



;i-28) M[9, X }=[ M ° ° r ) + ' ma ^ 



1 .An V 



where M is the differential operator 

(1.29) M [y] = -— div ( 6 2 k(6, X ) V^- 

Cy \ Cy 

with homogeneous Neumann boundary condition, and rriij are scalars given by the formulas 

(1.30) mil = J^Wf, 



e 

H(0) Cy 



1.31) m 12 = -— F7 ^—D x £, 



(1.32) m 22 



Note that m\ 2 = ^11^22 and mn > 0,m 2 2 > 0; hence, M^,^] is positive semidefinite. 
Furthermore, we have 



1.33) 



D e £\ ( Cy \ ( Cy \ ( D e S\ ( Cy/9 



D X SJ \D X E) \le + D x B)> \d x S J \& 



We easily check that ( I1.26P holds, and ( I1.27P has the form 



1.34) 



t \ _ /div (/c(6>,x)V6»)/cy\ (m n m i2 \ (c v /6 



XtJ V J \m 12 m 22 



ds 
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In component form, we have 



1.35) 9, = ^{&v(k(e,x)V0) + -^D x SD x A 



(1.36) Xt = "^V- 

To see that fll.35ft - fll.36ft coincides with fll.2ip - fll.22p . it suffices to take into account the 
formula 

c) c) 1 

^e = c v 9 t + — e Xt = div (*(0, x) V0) + T^y^x^ iV 7 = div (fc(0, *) V0) - Xt D X B . 

This shows that the model is compatible both with the standard principles of thermodynam- 
ics and the generalized thermodynamic formaliNotice thatsm introduced in [16]. Note that 
Mielke [22] recently elaborated the GENERIC approach for the phenomenology of thermoe- 
lastic dissipative materials. 

We prove an existence result for a suitable variational formulation of system fll.ll) - 
(11.41) . Using a Moser technique, we also show that the temperature variable 9 is globally 
bounded from above. The uniqueness result holds true for particular classes of potentials tp 
provided that the heat conductivity k in fjl.lj) does not depend on x, 7 = in (j 1.3ft . and 
cy and fi satisfy a suitable growth condition around 0. 

The paper is organized as follows. In Section [2J we state our assumptions on the data 
and our main results; in particular, global existence for a suitable variational formulation of 
fll.lft - fll.4p . In Section [3], we prove some auxiliary results related to the Lipschitz continuity 
of solution operators to differential inclusions. The proof will be developed as follows: the 
problem is approximated by partial time discretization, regularization and a cut-off procedure 
(cf. Subsection 14. ip . Suitable a priori estimates (cf. Subsection 14.21) allow us to pass to 
the limit with respect to the time step and regularization parameters, while the cut-off is 
removed by proving an upper bound on the absolute temperature (which is independent of 
the truncation parameter) by means of Moser techniques (cf. Subsection 14. 3p . The uniqueness 
result is proved in Section 



2 Main results 

In this section, we state our main results on solvability conditions for the system fll.lft - fll.4p . 
We start by introducing a suitable variational formulation; to this end, we consider a bounded 
and Lipschitz domain Q C M. N , N > 1 , and for t G (0, 00] we denote by Q t = f2 x (0, t) the 
open space-time cylinder and by E t its lateral boundary dQ x (0,t). We use, for the sake 
of simplicity, the same symbol H for both L 2 (fl) and L 2 (f2; WL N ) , while for an arbitrary 
integer d, H denotes the space L 2 (f2;IR d ). H and H are both endowed with the standard 
scalar product which we denote by (•, •) . The symbol V stands for the space i/ 1 (f2) , and V 
for its dual space, while the symbol V denotes the space if x (f2; M d ) , (•, •) being the duality 
V — V and V — V. Then, the following dense and continuous embeddings, where we 
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identify H (and H) with its dual space H' (and H'), hold true: V ^ H = H' ^ V , and 
V <^-> H = H' <^-> V . Finally, we rewrite the system (ll.lj) — (jl.4j) in the following variational 
formulation: 

(2.1) (dt(e(6,x)),z)+ [ k(6, x )V6-Vzdx+ [ 1 (9-9 r )zdA 

Jn Jan 

= - (X'(x)d t X + P (<p(x)) t + b[x}Xt) zdx Vze V, a.e. in (0, oo) , 
Jn 

(2.2) fi(9) Xt + A'(x) + ea'{x) + (P + B)dip{ X ) + b[ X ] 

+ e x (0,x) - ds x (6,x) 3 a.e. inQoo, 

where (12.21) has to be understood as an inclusion in R d with b[x] defined by (11.131) . and e 
and s are defined in fll.16p - fll.17l) . Letting (cf. (11.41) ) u := e(9 , Xo) ? we prescribe the initial 
conditions 

(2.3) e(0,x)(O) = u o , x(0) = Xo a. e. in Q , 
and suppose that the data fulfil the following assumptions. 

Hypothesis 2.1. (Existence) Let us fix positive constants C a , C\, k , k±, c, c, c\, (3, C , 
and assume that 

(i) ip : R d ->■ MU {+00} is a proper, convex, and lower semicontinuous function, T>(ip) is 
its domain; 

(ii) a, A G W^°°(2%)), |o»| < |A'(r)| < C A for all r G V(<p); 

(iii) k G W 1,00 (n x fi), = a.e. in Q x ft, G G ^ 2 '°°(X%) - 2%)), 

= G(-z) for all z G (V(ip) - V(ip)) ; 

(iv) k : R x P(<£>) — ► (0 + 00) is a locally Lipschitz continuous function such that < 
k < k(v, w) < k\ for all v G R and w G T>(<p) ; 

(v) The function \i maps [0, 00) in (0, 00) and the function 9 1— > j s bounded and 
Lipschitz continuous on [0, 00) with Lipschitz constant L^; 

(vi) cy : [0, +00) x V{ip) — > [0, +00) is a continuous function satisfying 

(2.4) Cy (0, X )=0, 0<c y (#,x)<c V #G(0,+oo), V X G%)i 

(2.5) c<c y (£,x) V(^x)G[l,+oo)x%); 

£y(0 x) 

(2.6) the function 9 \-t — is integrahle in (0, 1) for aii x e ^(^)- 

Moreover, for aii (0, x) G [0, +00) x £>(<^) there exists the gradient (cy) x (#, x) , and it 
hoids 

(2.7) \(c v ) x (9,x)\ < c lCv (9, X ), \(c v ) x (9,xi) ~ (c v ) x (9, X 2)\ < c x \xi - X2I 
for aii 9 G [0, +00), X , Xi, X2 e 2%). 
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Let e and s be defined by formulas ( jI.I6]) - (II.17]) and suppose that 

(2.8) 0< S (1, X )<C! V X 6%); 

(2.9) X i) - s x (9 2 , X 2)\ < ci (|0i - 2 1 + |xi - X2I) 
for all 61,62 G [0, +00), Xi, X2 G T>(<p). 

(vii) xo G V n L°°(fi) d . Moreover, for any C > set 

Z>c(vO = {X6 2%) : ^ £ € <Mx) : 161 < C}, 
and assume that Xo(%) G ^Coi^) a. e. in O; 

(viii) ^o, u G L°°(f2) fulfil u = e(0 o ,Xo) an d ^ofa) > a. e. in Q; 

(ix) 7 G L°°(dQ) is a nonnegative function; 

(x) 6 T G ^(Soo) is such that 6 r (x,t) > a.e. and log(0 r ) G L^E^); 

(xi) Let C G (W^'*(0, oo)) d be the solution to the differential inclusion 

(2.10) a(t)Ct + dcp(C)3g(t) a.e., 

with the initial condition ((0) = ( , ( G M d , and given data g G (L°°(0, oo)) d and 
o G Lj^ c (0, 00) such that < oo < ot(t) a.e. We assume that there exists a positive 
constant D > such that for all C > such that \g(t)\ < C , and Co G T>c(<fi) , we 
have 

(2.11) |G/-aC*)(*)l a.e.. 

We are now in position to state the existence theorem. 

Theorem 2.2. (Existence) Let Hypothesis 12.11 hold. Then there exists at least one pair 
(0, x) that solves system (I2.ip - fl2.3l) and such that 

(2.12) G L°°(g oo ) n L 2 (0, 00; V), (e(0, X ))t G ti^O, 00; V) , 

(2.13) 6(x,t)>Q a.e. inQoo, 

(2.14) xeTO^n^oojV), xt g L°°(c7oo) d ; 

3 C > : x(^, t) G £>cM a. e. in Qoo • 

Moreover, there exists a positive constant independent of t such that the following uniform 
upper bound hold: 

(2.15) 9(x,t) <0 for a.e. (x,t) G Qoo- 

Hypothesis 2.3. (Uniqueness) Assume that Hypothesis 12. II is satisfied and suppose more- 
over that 

(i) fc(0, x) = fc(0) for all G E and x € R d ; 
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(ii) Fix T G (0, oo) and suppose that there exists a positive constant R, depending only 
on C , a , and T such that the solutions d, ( 2 £ W 1,oo (0,T) to (I2.10p associated with 
data Coi,Co2 G ^cif) , &i,&2 G L°°(0,T), and with gi,Q2 G L°°(0,T) complying with the 
constraint 

(2.16) |&(t)|<C, i = l, 2, a.e. in(0,T), 

satisfy for all t G (0, T) the inequality 
(2.17) 



jf IC1-C2IO-) dr+|Ci-C2|(*)<i*(|Coi 



'+/( 


1 


1 






«2 



(t) + -^2|(r) dr 



(iii) Define c(6) := min{cv(?;,x) : X £ ^(v 9 ); u > ^} ai >d assume that 

Jo 

+00; 

(iv) The function v 1— >■ v 2 /u.(v) is nondecreasing in (0, +00) ; 

(v) There exists 9* > such that ^o( a; ) > f° r a - e - £ G f2, and 7 = on <9f2; 

(vi) Assume that c(9) > for every 9 G (0, 00) . 



c(v)fi(v) dv 



V 2 



Remark 2.4. Hypothesis 12 . II allows for physically meaningful choices of cy . We can choose, 
for example, a (d+1) -component model with K = {xi > 0, ^=1 Xi < 1}; Xo = 1 — Si=i Xi> 
</? = /a' (indicator function of the set K), and cy(#,x) = Si=i c *(^)Xij where q(#) behave 

asymptotically at and 00 like - — ^ a . Further examples of potentials ip complying with 
Hypothesis I2.1( xi) and Hypothesis 12.3( h) will be given in the following Section [3J 

We state then our last result regarding uniqueness and continuous data dependence 

for 

Theorem 2.5. (Uniqueness) Suppose that Hypothesis 12.31 is satisfied. Let T G (0, 00) be 
fixed. Then, there exists a positive constant 9(T) such that 

(2.18) 9(x, t) > 0{T) for a.e. (x, t) G Q T ■ 

Moreover, if (9i,Xi), (#2, X2) are two solutions to (I2.ip - fl2.3p in the sense of Theorem 12.21 
associated with initial data #01, X01 an d #02, X02 , respectively, and 9 = 9\ — 9 2 , x — Xi ~ X2 , 
Xo — X01 ~ X02 , 9 = Oqi — 9q2 , then, there exists a constant Ct > such that 



(2.19) / [ \9(x,t)\ 2 dx dt+ max [ \x(x,t)\ 2 dx<C T (\9 \ 2 H + \xo 

Jo Jn ^[o,T]J n V 



Finally, beside Hypothesis I2.3[ assume that 

(2.20) 9 eV. 

Then, the ^-component of the solution (9, x) to (I2.1]) - (j2.3]) has the further regularity 

(2.21) 9 G L°°(0, T; V), 9 t G L 2 (0, T; H) . 



11 



3 A differential inclusion 

This section is devoted to the description of some properties of solutions to general differential 
inclusions of the form ( I2.10p , which are used in the proof of Theorems 12.21 12.51 

First we provide some examples of functions f satisfying Hypothesis I2.1( xi). and 
we prove some further properties for space and time dependent differential inclusions that 
follow exactly from this assumption. Finally, we give examples of functions f satisfying 
Hypothesis I2.3l (ii). 

Examples of functions complying with Hypothesis I2.11 (xi). 

Proposition 3.1. The function if introduced in Hypothesis \2.1]f i) also satisfies Hypothe- 
sis \2.lV xi) in each of the following cases: 

(a) if d — 1 ; 

(b) if ip is the indicator function Ix associated with a closed, and convex set K cR d . In 
this case one has D = 1 ; 

(c) if <f(x) = f(M K (x)), where f : [0, / ) — > [0, +oo) is an increasing and convex C 1 
function such that /(0) = /'(0) = 0, fo > 0, and Mk is the Minkowski functional of 
K , a closed, convex set in M. d such that B r (0) C K C -Br(O) , defined by the formula 
M K (x) = inf {s > ; ±x e K} , for x £ R d . Then D = R/r . 

The proof of point (a) follows directly from [181 Prop. 3.4], (b) is obvious. Let us 
prove the point (c). In order to do that, we first need to prove the following auxiliary result. 

Proposition 3.2. Let f(x) be as in Proposition 13. 1 ]f c) . Then for every C > there exists 
C\ > such that for every x £ T>(f) we have the following implications 



We clearly have df(0) = {0}. For x ^ 0, take 7 £ (0, 1). Then, for rj £ df(x) and for all 
y £ T>(f) , we have 



(3.1) 
(3.2) 



<p(x) < Ci sup{|r/| : r] £ df(x)} < (R/r)C, 
(f(x) > Ci ^> inf { ] 77 ] : 77 £ df(x)} > C. 



Proof. We first prove the following equivalence 



rj £ df(x) (r) = cw, wedM K (x), c = f\M K (x))) . 



(77, x - (x - j(x - y))) > <f(x) - f(x - j(x - y)) , 



hence 



(t?, x - y) > - (f(M K (y)) - f(M K (x)) - j(M K (x) - M K (y))) . 
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77 

Letting 7 tend to 0, we obtain that — — — G dM K (x). Conversely, for w G dM K (x) , 

f'(M K (x)) K ' K 1 

we have 

(f(M K (x))w, x-y)> f'(M K (x))(M K (x) - M K (y)) > f(M K (x)) - f(M K (y)), 

which we wanted to prove. 

Let now C be a given positive constant. For all w G 8Mk{x) and 1 ^ we have 
(1/R) < \w\ < (1/r). From this we deduce that if f'(M K (x))\w\ < C, then f'(M K (x)) < 
CR, and so M K (x) < (D'^CR) and <p(x) < /((/' • We can choose d = 
/((/') _1 (C-R)) , and (13.21) is proved. Suppose now that f(M K (x)) < C±. Then we have 
M K {x) < {f)- l (CR), hence f'(M K (x)) < CR and f(M K (x))\w\ < (R/r)C , and (f37TJ) is 
proved. ■ 

We conclude the proof of Proposition 13.1( c) by proving the following Proposition 13.31 

Proposition 3.3. Let if be as in Proposition \3.V( c). Then Hypothesis \2.W xi) is satisfied 
with D = R/r. 

Proof. Consider some ( satisfying inclusion (I2.10p with initial datum ( . Then, the 
following equality holds true for all t G (0, 00): 

Ht)( t \ 2 + \g(t) - a(t)Ct\ 2 + 2a(t)<p(() t = \g(t)\ 2 ; 

hence, we immediately deduce that, for all t G (0, 00), 

(3-3) y{Q t = ^ {\g(t)\ 2 - \g(t) - a(t)( t \ 2 - \a(t)( t \ 2 ) . 

By assumption, we have \g{t)\ < C and g(t) — a(t)( t G d(f((). In view of (13. 2p . there is a 
constant C\ such that 

Hence, 

(3.4) f(C)MC) ~ Ci) + < VtG(0,oo). 

Integrating from to t and using the assumption £ G U c (<p), it follows that </?(£) (t) < C\ 
for all t G (0, 00) . Then, using (13. ip . together with the fact that \g{t)\ < C for all t G (0, 00) , 
we finally obtain that 

(3.5) \g(t) - a(t)£ t \ < -C and \a{t)( t \ < ( - + l) C, 

r \ r J 

which concludes the proof. ■ 
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Properties of the solution mapping under Hypothesis l2.lT xi). 

Proposition 3.4. Let us consider the solutions Ci, C2 G W 1,o °(0, 00) to (I2.10p associated 
with data £ i, C02 G ^clv 9 ) , «i, «2 G -kto c (0> 00 ) > aDC ^ with g\,g% G L°°(0, 00) complying with 
the constraint 

\9i{t)\ <C, i = 1, 2, a. a, 

and let Hypothesis l2.II( xi) hold. Then there exists a constant L such that for every t G (0, 00) 
we have 

(3.6) |Ci " ClW<ICoi- Coal + ^jT( 



1 



1 

Q!2 



Ifi 1 ! - #2 



dr 



Proof. Test the difference of the two inclusions (12.101) by £i — £ 2 and divide the resulting 
inequality by aii . Then we obtain for a.e. t G (0, 00) : 



<Ci-C2,Ci-C 2 > < 



1 

Oil 



1 

«2 



("2(2,0 - (2)! + \ (gi - 92, Ci - Ca) 





1 


1 








+ 1^1 ~02^ 




«1 


OL 2 





Using the bound for | CK2C2 1 (cf. Hypothesis I2.1l (xi). (12. lip ), we get 

from which (I3.6P immediately follows by integrating over (0,t). I 

Proposition 3.5. Let Hypothesis \2.lV xi) hold, and let ( n and ( be the solutions of f)2.10p 
corresponding to the data (g n , a n , ( 0n ) and (g, a, ( ), respectively, with \g n (t)\ < C , 
Con G T>c(f)- If {Con} converges to ( in R d , {g n } converges to g and {a n } converges to 
a in L 2 (0,T) for some T > 0, then {( n } converges strongly to ( in L 2 (0,T). 



Proof. Test (I2.10p . written for ( n , by ( n in order to obtain 



9n ~ «nCn ) Cn = "^(Cn) • 



Now set r] n 



9n 



'OLr 



2-v^nCn- Then, by straightforward computations, we obtain that 



1 OLr, 



\Vn\ 2 = 4 ^^(Cn) • 



We know, by Proposition I3.4[ that converges uniformly to £ and that (p is Lipschitz 

continuous on T>o(<fi) ■ Hence, integrating over (0,t), we obtain that |?7nU 2 (o,T) — > \v\l 2 {o,t)- 
Since we know that n n — » n weakly in L 2 (0, T) (since ( n — > £ weakly in L 2 (0, T) ), we infer 
rj n — > T) strongly in L 2 (0, T) , which is sufficient in order to conclude the desired convergence. 
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Examples of functions complying with Hypothesis l2.3T h). 



Proposition 3.6. The function ip introduced in Hypothesis \2.V( i) satisfies Hypothesis \2. 3( ii) 
in each of the following cases: 



(a) if d = 1 ; 

(b) if, for any C > , ip is 

(c) if ip = Ik, where K is either a polyhedron or a smooth 
interior. 



a C 1 -function with Lipschitz continuous derivative on T>c(<p); 

convex set with nonempty 



Proof. The proofs of (a) and (c) follow respectively from [T8J Prop. 3.4] and [3, Thm. 7.1, 
p. 88]. We briefly show here how to proceed to prove case (b). Let us consider solutions 
Ci, C2 G W^iO^T) to (12301) associated with the data Coi,Co2 e V c (ip), a ± ,a 2 e L°°(0,T), 
and with g±,g 2 € L°°(0,T) complying with the constraint 

\9i{t)\<C, i = l,2, a.e. in(0,T). 



By Hypothesis I2.1f xi). (i,( 2 remain in T>r>c(f)- Using the Lipschitz continuity of ip' on 
U DC (<p) , we obtain that there exists a positive constant Q such that the following inequality 
holds true a.e.: 

«i|Ci - C2I < |«i - CK2IIC2I + Q\Ci - C2I + Ifl-i - 92]- 

Dividing by ot\ , and using the bound for |ck 2 C2 I , from Hypothesis 12 . 1 ( xi) (cf. (12. lip ), we get 



(3.7) 



ICi - C2I < (C + DC) 



1 
1 

a 2 



1 

a 2 



Q\(i -C2I + \gi-92\ 
ICi - C2I + \gi-92\ 



for some positive constant M (depending on C,D,Q). We can rewrite this inequality in 
the following convenient form, for t G (0, T) , 



d 
dt 



Mt f 


1 


1 




«1 


«2 



+ \9i ~ 92 



(e- Mi |Ci-C2|)<e- Aft ( 
Integrating over (0,t), and using the previous inequality (13. 7J , we deduce 

ICi — C2K*) < ^ I IC01 — Coal 



1 


1 






1 


1 








+ \gi - 92] -f 


■a 






+ \gi 




a 2 








a 2 





dr 



for some positive constant L (depending on C,D,Q). Integrating once more in time we 
arrive at the desired inequality (I2.17p . ■ 

Proposition 3.7. Let f : [0, fo) — > [0, +00) be an increasing, convex function with locally 
Lipschitz continuous derivative, f(0) = f'(0) = 0, and let K be a closed, convex set of class 
C 1 ' 1 such that B r (0) C K C Br(0) . Then, <p(x) = /(M^(a;)) has a Lipschitz continuous 
derivative on T>c{<p>) for any C > 0, i.e., property (b) in Proposition 13.61 is satisfied. 
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Proof. Let C > be given. We denote V c f = {s e (0, /„) : f(s) < RC} , and let 
Lc be the Lipschitz constant of /' on T>cf ■ For x G T>ci}P) we have |y?'( x )l — @ > hence 
f'(M K (x)) < RC , that is, M K (x) e X>c/- We now estimate the difference |</?'(x) — <p'(y)\ 
on T>c(<p) ■ Assume first that x ^ 0, y = 0. Then 

- = \<f/(x)\ = f{M K {x))\M' K {x)\ < \f{M K {x)) 

< -L c M K (x) < \l c \x\ = -%-\x - y\ . 

Consider now the case x ^ 0, y ^ and set Jk(x) = M K (x)M' K (x) , Jk{u) = M K (y)M' K (y) . 
The mapping is Lipschitz continuous on M d (with Lipschitz constant Lj ) (see p], Section 
5.2]), and we have 

y{x) - <f/{y)\ = \f{M K {x))M' K {x) - f{M K {y))M' K {y)\ < f '[f K ( ( f \J K (x) - J K {y)\ 

M K (x) 

+ IM^^I ^^V k^) - M K (y)\ + \M' K (y)\\f(M K (x)) - f'(M K (y))\ 
< L c (^Lj + \x-y\, 

from which the assertion follows. ■ 
A relevant case for applications is, for example, <p(x) = — log(l — Mj^x)) , see [IT] . 



4 Existence of solutions 

This section is devoted to the proof of the existence result stated in Section El We use a 
technique based on approximations, a priori estimates, and passage to the limit. 

Let us first write down our equations (12 .ip and (12.21) as 

(4.1) ((e(9,x)) t ,z)+ [ k(9, X )Ve-Vzdx+ [ 1 (9-9 r )zdA 

Jn Jan 

= - I (^(x)Xt + H t P(x))t + b \x]Xt) zdx VzeV, a.e. in (0,oo), 
Jn 

(4.2) n(0) X t + (P + e)Mx) i -A'(x) - Mix) - M - <hc{e, x) + 0s x (9, x) 

a.e. in Q^. 

4.1 Approximation 

Assuming Hypothesis 12.11 to hold, we proceed as follows: first we extend the domain of 
definition of cv(9,x) by putting cv(9,x) = cv(\9\,x) f° r X) £R X ^(v 9 ); an d set 

e(0, X )= [ c v (C,x)dC for(^ lX )6Mx%). 
Jo 
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U{\6\) x\e\<Q 
\v(e)(\9\- e ) if |d| > g 



We now fix a truncation parameter g > 1 , which will be determined below, and define 

Jl e (9) 

We fix an arbitrary T > 0, and split the interval [0,T] into an equidistant partition = 
to, ti, . . . , t n , tj = jT/n for j = 0, 1, . . . , n, n G N, with the intention to let n tend to oo . 
We choose sequences {9o tn } n G V and {#r>}n £ H^ 1,2 (0, T; L 2 (9f2)) of approximate data 
such that #0,71(20 — V n a - e -5 #r,n(^)*) > V 77, a - e -j #o,n — ► #0 strongly in H, and #r> - ► 9 r 
strongly in L 2 (0, T; L 2 (dtt)) . 

An approximate solution (9 n , Xn) will be constructed successively in intervals [tj-i, tj] 
for j = 1, . . . , n. Assuming that it is already known on [0, , we define 



(4.3) Xn(x,t) = Xn(x,tj_i) xett, tefa-utj), j = l,... 

Xn forte[0,ti) 
n (x,t) = < n rt s -i 



,n, 



7 1 



/ n (z, r) dr for t G fo-i, tj), j > 2 . 

Jtj-2 



With this notation, we then state the following approximating problem. We use only the 
index n for the variables here (omitting the g dependence), for simplicity. 

Problem (P) m . Find two functions 6 n e H l {0, T\ H) n L°°(0, T; V) and x« e L°°(fi x 
(0,T)) d , 9tXn G L°°(fi x (0,T)) d , such that 9 n > e n a.e. in Q T for some e n > 0, Xn G 
T>c((fi) , and for all £ G (0, T) and z G V, we have 

(4.4) 

jf <9 t (^n(t) + e(0 n (*), Xn(0)J ZdX + j Q k ^n(t), X n (t))V9 n (t) ■ Vz dx 



+ / 7 (0 n (*)-0r, n (*))*<L4 = - / ((A'(xn)(t)+MXn](t))ftxn(t) + /3^(^( X n(t))))^dx 
Jan 

(4.5) 

Ji e (9 n (t))d tX n(t) + ((3+ \9 n (t)\)MXn)(t) 3 -\'(Xn)(t) ~ \9 n (t) \a'( X n) (*) - &[**](*) 
-e x (^(t),X„.(t)) + |0n(*)l4(^n(*),X«(<)) a - e - in 

with initial conditions 



(4.6) e n (o) = e , n , xn(o) = Xo . 

Lemma 4.1. Under Hypothesis \2.1\ for each g > and n G N tie Problem (P)( n ,e) ^ as 
a unique solution (9 n , Xn) with the required properties. 
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Proof. On each interval (t,_i, tj) , we can proceed as in the proof of [T9J Thm. 2.2, p. 290]. 
We test a Galerkin approximation of ( 14.41) by the approximation of d t 9 n . The estimates are 
sufficient to pass to the limit in the Galerkin scheme and to obtain a solution on each interval 
(tj-\,tj). We only have to check that the initial conditions at ti,t2, ■ ■ ■ are well defined. 
Indeed, since on each interval (tj-x,tj) we have 9 n £ H x (tj-\,tj\H) n L°°(i,_i, tj] V) , we 
also obtain that t i— > 9 n (t,-) is weakly continuous in (tj_i,tj) for every j with values 
in V. Moreover, \n is strongly continuous with values in L°°(Q) d , and there exists a 
positive constant C (independent of n) such that Xn(t,-) £ Vc{<f) on (tj-i,tj) for every 
j = 1, . . . , n. Hence, we can define the initial conditions at t = tj by 9 n (tj) = 9 n (tj—) . ■ 



4.2 A priori estimates 

In this subsection, we perform suitable a priori estimates (independent of n) for the solution. 
In the following, we will denote by C any positive constant that depends only on the data 
of the problem but may vary from line to line. In particular, it will not depend on the 
truncation parameter g and discretization parameter n. If such a dependence takes place, 
we use the symbol C e for a constant that depends on g, but not on n. Again, the same 
symbols will denote constants that may differ from line to line. 

Let us, for simplicity, in this subsection occasionally omit the indices n and write 
simply 9,x instead of 9 n ,Xn if no confusion arises. We denote (note that 9 n > and so 

e = e) 

(4.7) u n {t) := -9 n (t) + e{9 n {t), Xn (t)) for t £ (0, T) . 



Estimate for Xt • Equation (14. 5p is of the form f)2.10p with 

(4.9) 9 (f) = <[», X \ := -j^g (6°'(X) + Ax) + »M + ^(9, x) - »««(», x)) ■ 
First, let us note that, using (12.41) . (12. 7p . and ( 12. 8p . we get 

(4.10) \) < /'■"" ,s ' t, K-M^ d? < Cl £szMl i( + Cl ^' | d? 

< cis(l, x) + ciclog £> < c^ + ^clog^. 



Hence, owing to Hypothesis 12. 1[ we have 



x]\<C a + -AC x + C b ) + -^- sup |( Cv ) x (e )X )| + -^--|4^ x )| 

P P + v 0<£<9 P + V 

< C a + - (C x + C 6 ) + ci sup |c y (£, x)| + cl + ciclog g 

P 0<£<6» 

< C CT + -g (Ca + C b ) + cic + c\ + ciclog g , 
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where C& denotes here the upper bound for the operator b defined in (j!.13p . Let us set 

(4.11) Ce, e := C a + i (C x + C b ) + c x c + c\ + ciclog e . 

Then the conditions of Hypothesis I2.1f xi) are satisfied with the choice C = ma,x{Ce tQ , Co} , 
where Co is defined in Hypothesis I2.1( vii) and Ct )B is defined in (14. lip . Hence, we obtain 
the following estimates on x'- 

(4.12) \Xn\L°o(Q T )+\d t Xn\L~(Q T ) + \dt((p(Xn))\L™(Q T ) < C(l+log^) 2 , 

where now C is a constant independent of g. 
Estimate for 9. Taking z = 6 n in (I4.4p . we get 

(4.13) / dt (-8(t) + e(9(t), X (t))) 6(t) dx + [ k(d n (t), *„(*)) V0(t) ■ V0(t) dx 



j(6(t)-6 r , n (t))e(t) dA 
an 

- ! (X'(x)(t)d t x(t) + (3d t ( V ( X (t))) + b[x}(t)dtX(t))0(t)dx. 
Jn 

U(9,x)= f c v (v, X )vdv for(0,x)e(O,oo)x2%). 



o Jn 



Define 



We have d t U{6,x) = 0d t e(6,x) + (d x U - 6d x e)xt- We integrate (T4~T3D from to t and 
rewrite the first term as follows: 

\(h(T) + e(9(r),x(r))\ 6(r)drdx 

[^(t) + U(6(t), X (t)) S ) dx-^^ 2 (O) + f/(0(O), X (O))) dx 

1 [ (d x U(e(r),x(r))-e(T)d x e(e(T),x(r)))xt(r)dTdx. 
o Jn 

There exist two constants C\,Ci such that U(6,x) > C\0 2 — C%. Hence, by (14.121) and the 
Gronwall's lemma, we obtain 

(4.14) |IM£»(o,T i V)nL<»(o J T;£0 ^ C ( X + lo S • 
By comparison, we also deduce that 

(4.15) \\d t u n \\ LH o,T;V>) <C(l + \ogg) 2 . 
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Estimate for Vx- The function 



(3 + 9 



is Lipschitz continuous in R due to Hypothesis I2.1f v) and, with the help of the mean value 
theorem, it is straightforward to deduce that 

(4.16) X i] ~ Ifo, X2]\ < k'(Xi) - cr'(X2)\ + y l^i " Oal + 4|A'(Xi) ~ A'(x 2 )| 



/3' 



+ §rl*i - fc| + ^[Xi] - &NI + pl^i - 02 1 + 



/? 2 



/3 



^/ lci|Xl - X2 ' de 



4(^x2) I 



+ -|^(0 1 , Xl )-4(02,X 2 )| + 



i^ 2 



-1^1 - ^2 



< k'(xi) - ^'(X 2 )l + -g (|A'(xi) - A'(x 2 )| + \b\xi] - b[ X 2}\) 

+ j 2 (c„p + c A + c b ) \e x - e 2 \ + j\xi- x 2 | + - (cic + Cl c) ^ - 2 | 
+ 1 - 2 | + Ixi - X2I) + ^ (c? + Cl log q) \e x -e 2 \. 

Hence, we can apply estimate (13.61) in Proposition 13.41 to ( 14. 5 p with (for 



(3 + 9 

(Ml) 



(y,t), 



ai(t) 
a 2 (t) 
9i(t) ■ 



g 2 {t)=£[9(y,t), X (y,t)] 



\(3 + 9j 

M),xOM)] 



1 



0</(x) + A'(x) + &[x]+ex(0,x) 



(3 + 9 

-9s^(0, X ))(x,t) 
1 



Mix) + A'(x) + b[ X ] + e x (9, x) 



and 

Hence, we obtain 



(3 + 9 
-6s^e,x))(y,t), 

Ci = X(x, t), C2 = X(V, t), x, y G ft. 



(4.17) \ X (x,t) - X (y,t)\ < \Xo(x)-Xo(y)\ + L e { J \9(x,s)-9(y,s)\ ds 

+ / (\ b [x](x,s) -b[x](y,s)\ + \x(x,s) - x(y,s)\) ds 
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where L e depends on L, L M , and the constants on the right-hand side of (14.161) . Now, 
recalling (jl . 13j) . we have, by ( I4.16p . that 

(4.18) \x(x,t)- X (y,t)\<\ Xo (x)-Xo(y)\ + L g [ \e(x,s)-6(y,s)\ds 



o 



+ 2L e / \K(x,z)(G'(x(x,s)-x(z,s))-G'(x(y,s)-x(z,s)))\ dzds 
Jo Jn 

+ 2L e / \G'(x(y,s)-x(z,s))(K(x,z)-K(y,z))\ dz ds 
Jo Jn 

+ L g / \x{x,s) - x(y,s)\ ds. 
Jo 

Thus, in view of Hypothesis I2.1l (iii). we obtain that 

(4.19) \x(x,t)-x(y,t)\<\xo(x)-Xo(y)\ + L Q [ \9(x,s)-9(y,s)\ds 

Jo 

+ L e {2L b + l) \ \x(x,s) - x(y,s)\ ds + 2L g L b / \k(x, z) - n{y, z)\ dz , 
Jo Jn 

where L b is a constant depending on the Lipschitz constants of G and G' , ||re||z,°°(f}xn) > 
and T. From f)4.19p . using the assumptions xo £ V and k G W 1,00 (Q x Q) (cf. 
Hypothesis I2.1( iii)). we immediately deduce that 

\V X (;t)\<C e (l + jT (\V6(-,s)\ + \V X (-,s)\) ds^j a.e. in Q . 
Now, with the help of Gronwall's lemma, we infer that 

(4.20) |VxM| <C e (l + J |V0(.,5)|ds^ a.e.infi. 
Using finally (14.201) with (14.141) . we get the desired estimate 

(4-21) ||Xn|U°°(0,T;V) < C Q , 

where C e denotes a positive constant depending increasingly on g. From the definition (14. 7p 
of u, it also follows that 

(4-22) ||«n|U 2 (0,T;V) < G Q . 

4.3 Lower and upper bounds on 9 

In this subsection, we first prove a bound for log^ entailing the strict positivity of the 
absolute temperature (in the limit when n — > oo). Then, we prove a (time dependent) 
upper bound holding true for the solution component 9 n for n fixed, which enables us to 
proceed with the Moser iteration procedure in order to prove a uniform (independent of time, 
of n, and of g) upper bound on 9. This permits us to remove the truncation parameter 
and to conclude the existence proof. Finally, we will prove a lower bound (independent of 
n ) on 9 holding true under the additional Hypothesis 12.31 that we will use for the proof of 
uniqueness of solutions. 
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Estimate on log 9. Let us rewrite equation (I4.4p . by using (14. 5p . in the following form, 
for all z G V , 

(4.23) / d t (-6 + e(6,x)) zdx+ [ k(9 n , x n )V0 • V^dx + f 7 (0 - 9 r , n ) zdA 

Jn \ n J Jn Jan 

TtMXt+OxtR(0,x)z dx, 



where 

R(6, X ) ■=cr / (x)-s^e,x)+^ eeW, \£(x,t)\<C&.e., 

C being defined in Hypothesis 12. ll fxi) . We prove now an estimate on log^ in L 2 (0,T;V) 
by taking in (14.41) z = T(9), where 

(4.24) T^.-^VJ'I 

V 6 J \Q for#>l 

Notice that we are allowed to perform this estimate, with fixed n, because 9 > e n > a.e. 
for all n (cf. Lemma FOi) . We get, using equation (14.51) and Hypothesis I2.1( iv). 

(4.25) A E (9, X )- J Xt^Mx&xmt) d^dx + k \\V(\og9)-\\ 2 H 

+ / (9-9 r , n )T(9)dA< f ]l e (9)\xt\ 2 T(9)dx+ [ 9T(9) Xt R(d,x) dx , 
Jan Jn Jn 

where 

E{d,x)= ~LL cv ^ x h i ~\) d ^°- 

Note that the first term on the right-hand side of (14. 25p is nonpositive, while the other term 
can be estimated using estimates (I4.10p . (I4.12p . and ( I4.14p on our solution (9, x) ■ Regarding 
the second term on the left-hand side in ( I4.25p . using (12. 7p in Hypothesis I2.1( vi). we obtain 
that 



< CiE(9,x)\\dtX\\L~>(nr- 
Moreover, we treat the boundary integral in the following way: 

(4.26) / (9-9 T , n )T(9)dA> [ t>{9) - V(9 r , n ) dA, 

Jan Jan 

where ^ is defined as 

f*-bg* for£<l 
\l for 9 > 1 ' 

and is a convex function on [0, +oo) such that ty'(0) = T{9) for all 9 e [0, +oo) . Using 
estimates (14.121) and (I4.14p . and Hypothesis 12.1 f x), we deduce that 

(4.27) E(9(x,t), X (x,t)) + k [ t \\V(\og9y\\ 2 H dt+ [ [ ^(9)dAd^ 

Jo Jo Jan 

<cjl+ [ E(9,x))+E(9 ,n,Xo)+ [ [ ^(^r,n)dAde 
\ Jo J Jo Jan 
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for a.e. (x,t) e Qt- Now, in view of (12.81) . we have that 



r#0,n 

(4.28) E(6 0>n , X o) = - I I c y (^Xo)(l-7) d Z dx 




n Ji 



< [ f c v ^,Xo)(l-l] d£dx< [ s(l, Xo ) < c x |n| . 
Jo Vs / Jn 

Moreover, we can estimate the term J Q * J dn \P(#r,n) dA dr using Hypothesis 12. l( x) as follows: 

(4.29) / / *(0r,n)dAde<C(||0 ri n|U-(E t ) + ll(log(^))"|Ui( Et )) < C - 
Jo Jan 

Using a standard Gronwall's lemma in (14.271) . together with the estimates (I4.29P and (I4.28p . 
we obtain the desired bound 

|| (log n )~ 11x2(0^) < C e , 
which, together with estimate (14. 14j) for 9 in L 2 (0,T; V) , gives 

(4.30) \\log9 n \\ L 2 {0>T . v) <C g . 

Lower bound for 9 under Hypothesis 12.31 (iii)— (vi). Let us consider relation ( I4.23p . 
Notice that, by the previous estimates, it follows that there exists a positive constant R Q 
such that \R(9, x)\ < Rq a.e. in Qt ■ Hence, for every z G V such that z > a.e., we obtain 
the following inequality (notice that by virtue of Lemma I4TT| we have here 9t G L 2 (0, T; H) ): 

(4.31) jf [U t + c v (9, X )9?j z dx + jf fc(0 n , xn) V0 • dx 

We now compare this inequality with the following ODE: 

R 2 a w 2 

(4.32) c(w)w t = -t^V^' = w ° > 

where c is defined in Hypothesis (I2.3p (iii) and wq = mm xe n 9 (x) > (cf- Hypothe- 
sis ||v)). 

Notice that the solution iw is decreasing and does not vanish in finite time, due to 
Hypothesis I2.3( iii). The function w does not depend on x, hence we may add to the ODE 
in ( I4.32p the term — div (k(9 n , Xn)Vti?) , which is equal to 0. Using the fact that w t < and 
c(w) < - + cv(w,x), we obtain, subtracting (I4.3ip from (I4.32p . the inequality 

^ + c v {w,x) \ w t - Q + c v (9,x) \ S^j zdx 

f - R 2 f ( 9 2 w 2 \ 

+ / k(9 n ,Xn)V(w-9) ■ Vzdx < -f \ ttttx- z dx. 
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We now take as test function z = H £ (w — 6) , where H e is the regularization of the Heaviside 
function H, 

'0 ifw<0 
H e {v) = { v/e ifve (0,e) . 
1 if v > e 

By virtue of Hypothesis 12.31 (iv)-(v), we deduce 

n \ + c v {w, x)J w t - Q + X)J 0*J # e (™ " 0) dx < 0, 
and we can pass to the limit in this inequality for e \ , getting 

J((^ + Mw, x)j w t-(~ + MO, *)) 0t) # (™ - ^) dx < , 

that is, 
(4.34) 

Notice now that, by Hypothesis I2.1l (vi) (cf. ( 12. 4 p and (12. 7p ). we have 

\ e x( w 'X) ~ e x (0,x)| < max(c y (r,x))k - 0| < c|u> - 0\ . 

T<WQ 

Integrating (14.341) over (0,t), and using the boundedness of Xt m L°°(Q T ), we obtain, by 
the choice of the initial data 9 and u> , that 

(4.35) J (^(J^w + e(w,x)^ ~ (J^Q + e(6, X )^ (t)dx < c J {w - 6) + dx . 
For w > 8, we have 

'("■X) -«(»■*) = JTMr, X )dr > £cMdr > 1 r- = 
w — 6 w — 6 w — 6 w Jo 

The function C is nondecreasing, C(w) > for w > 0. Hence, 

(4.36) (^w + e(w,x)j - (^9 + e(6,x)j > C(w)(w-6) forw>9. 
Inequality (I4.35P then yields 

C(w(t)) [ (w-6) + (t)dx < C [ [ (w-6) + dxds 
Jn Jo Jn 

for every t G (0, T) . From Gronwall's lemma we conclude that 

(4.37) 6 n (x,t)>w(t) a.e.mQ T . 
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Upper bound for 9. Let us denote the right-hand side in (14.41) by 

M(9, X ) ■= ((A'(x)(t) + b[ X }(t))d tX (t) + Pdt(<p(x(t)))) 

which, due to the previous estimates is bounded by a positive constant, say, M e . Then, we 
compare the inequality, for all z e V , z > a.e., 

(4.38) J d t (U(t) + e{9{t), x(t))^j z dx + j k{9 n (t), *„(*)) V0(t) ■ Vz dx 

+ f 1 (9{t)-9 v . n (t))z&A< ! M e zdx, 
Jan Jn 

with the following ODE: 

(4.39) -v n = M e , v(0) = v , 

n 

where v = max{sup 9 , sup # r } (cf. Hypothesis 12.11 (viii) (x)). 
Then we have that 

v n (t) = v + M e nt . 

We proceed as above, adding to the ODE the term — div (k(9 n , Xn) Vf n ) , which is equal 
to 0, and subtracting (14.391) from (14.381) . We test the resulting inequality by H e {6 — v n ) 
(cf. pi| ). Using the fact that (1/rc + c v (w, x)) v n > and that 9 E H 1 ^, T; H) , we can 
let £ tend to , getting 

/J((^ +e M-6" +e H) +di 

(e x (e, x) - e x (v, x)) XtH(0 - v)dx < . 



Using now the Lipschitz continuity of e x (cf. Hypothesis I2.1f vi)) and the boundedness of 
Xt , we obtain that 



— v ) + dx 



/ fi l((^ + ^^)"fr +e( "' x) )) dx - c ' 

Integrating over (0,t) and using the choice of the initial condition v , we infer 

[ (9- v) + (t) dx < C e [ [(9- v) + dx dr, 
Jn Jo Jn 

and, applying Gronwall's lemma, we get the desired upper bound 
(4.40) 9 n (x,t) <v n (t) a.e.mQ T . 
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Moser estimate. In order to conclude the proof of existence of solutions to (I2.1l) - (l2.3p . 
it remains only to prove that the 9 component of the solution (9, x) is bounded from above 
independently of T, g and n. 

This will enable us to choose g sufficiently large in such a way that in this range of 
values of 9 we have s x = s x . To this end, we perform the following Moser estimate. 

We will make repeated use of the well-known interpolation inequality (cf. [5]) 
(4-41) \\v\\ H < A (v\\Vv\\ H + V~ N/2 h\\ L i (u) ) , 

which holds for every v G V and every 77 G (0, 1), with a positive constant A independent 
of v and rj. 

Following the ideas already exploited in [TjJJ Prop. 3.10, p. 296], for j G N, we choose 
in g3D z = ((0„ - e R ) + f 3 - 1 G L 2 (0, T; H l (tt)) , with 9 R = max{0 r - 1, 0, 1} , where 

\@r,n(x,t)\ < 0r a.e. in St, |^o,n(^)| < © a.e. in Q. 

We know that z G L 2 (0, T; H\n)) , due to the upper bound on 9 n proved in (14. 40 p . 

Here below, we denote by Cj, i = 1,2,... some positive constants that may depend 
on the data of the problem, but not on j, T, g, and n. We omit again the indices n in 9 n , 
Xn , for simplicity. 

Now set u = (9 — 9 R ) and take z = {u + ) 23 ~ l in (I4.4p to obtain that 

(4.42) (U t {t) + (e(9(t), X (t))h (^f" 1 ) + J k(0 n (t), Xn{t))Vd{t) ■ V^) 2 ^ 1 dx 

+ / 1 {9{t)-9 Ttn ){u + ) 23 - 1 dA 
Jan 

= ~ [ (\'(x)(t)Xt(t) + P(<p(x(t))t)-b[ X }(t)xt(t)) (u+f-'dx. 
Jn 

Our aim is to prove that there exists a positive constant C* (independent of T, g, and n) 
such that 

(4.43) \\9{t) || £oo(n) < C* (1 + log gf +2N for a.e. t G (0, T). 
The first term on the left-hand side can be rewritten as 

(4.44) /I^) + (e(£(t), x (^ 

\n I n 

+ (e x (0,x)Xt,(u + f- 1 ). 

Let us deal with the second term on the right-hand side in (14.441) . using Hypothe- 
sis EH(vi) (notice that by (12. 5p . we have cy > c in the set where 9 > 9 R > 1) as follows: 

(4.45) (^cv^xX^) 2 ^ 1 ) = f t E ^x) - J q ( K j\c v ) x ^ + 9 R ,x)(er~ 1 <%) Xtdx, 
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where 
(4.46) 

cTi [ (u+f dx<%(«,x) := / f 'cvit + e^xXef- 1 d£dx < c2~' [ (u + f dx. 
Jn Jn Jo Jn 

Then, using (14. 12p and Hypothesis I2.1( vi). we infer that 

(4.47) ^jf (cy) x (e + ^,x)(e + ) 2J_1 d£j X*da; < d(l + log g) 2 E 3 ( U) X ). 

Moreover, by Hypothesis 12. 1[ we obtain the inequality 

/ k(9 n , xn) VmV ((u+f- 1 ) dx + [ 7 (u - 6 r , n + R ) (u+f- 1 dA 
Jn v J Jan 



, 2* -1 

> k 



2V- 



2 



!> 



" dx+ I 7 ((u + f -(u+f- 1 ) dA. 

an 



Now, set $j = (m + ) 2J 1 . Regarding the terms on the right-hand side in (I4.42p and the last 
term in (I4.44j) . using (I4.12p and Hypothesis I2.1( vi). we realize that 

- f (X'(x)(t)Xt(t) + ^(^(x(t))) t -b[x}(t)xt(t)-e x (e, X )Xt) (u+f-'dx 
Jn 

<C 2 (l + logg) 2 (J\^\ 2 dx + l\ . 

Let us now set = Ej H / |$j| 2 dx. Then, with the help of Holder's and Young's 

n Jn 

inequalities, we deduce that 



dt 



E >>x) + k0{ ZJ ] [ |V%| 2 dx+ / 7 |^| 2 dA 
/ Jn Jan 

< (1 - 2~ j ) [ 7 dA + 2~ J f 7 dA + d(l + log gfE^u, X ) 
Jan Jan 

+ C 2 (l + \ogg) 2 / (|$,| 2 + l) dx. 
Jn 

Multiplying the above inequality by 2 J ' , in view of the upper bound for Ej(u, X ) in (14 .46 p . 
we find out that 

(4.48) 2 j -^-E«(u, X ) + 2k [ \V$ j \ 2 dx+ [ 7 |d> i | 2 dA 

at Jn Jan 

< yc 3 (i+\og Q ) 2 (i+ jj^\ 2 d x y 

We now use the interpolation inequality (I4.4ip and note that, thanks to estimate ( I4.14p . we 
have 



||$ ii 



HI^(O) = (j n U+dX ) <C4(l + lcg£) 4 , H^llinoj = ||%-1|| H 
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Thus, we derive the inequalities 

(4.49) I \^\ 2 dx< 2A 2 (r ] 2 \\V$ 1 \\ 2 H + r ] - N C i (l + \ogQ) 4 ) , 
Jn 

(4.50) / |$/ dx< 2A 2 (^WV^fH + V~ N \\^3-i\\h) for J> 1 - 
Jn 

For j = 1 , we infer from (I4.48[) and (14.491) that 

2^r(«,x) + 2A;o / |V$ 1 | 2 dx+ / i\^\ 2 dA 
at Jn Jan 

< C 5 (l + \ogg) 2 (l + r) 2 J |V$i| 2 dx + r]- N (l + \ogg)^ . 
Choosing 77 = v^o/ {\JC^{1 + log g) 2 ) , we find that 

(4.51) 2-%-E?{u, X ) + ko f \V$ 1 \ 2 dx+ [ 7 |$i| 2 dA<C , 6 (l + logf ? ) 

at Jn Jan 

For j > 1 , we get 

< 2^CV(l + log£ ? ) 2 ri + ^jTlV^dx + ^llVillir)- 



6+7V 



Choosing 77 = y/k^/(y/2^(C 7 (l + log £>) 2 )) , we conclude from ffl~4"gj) and fl4"3U|) that 

(4.52) 2»"^£y(u l x) + MV* i ||| r + / 7^f dA 

dr Jan 

< 2^ +1 )c 8 (l + log£) 2 (1 + (1 +log^) Ar ||$ J _ 1 ||y . 



By assumption, we have ||$j(0)||# = 0. Hence, integrating (I4.5ip and (I4.52p with respect 
to time and using the lower bound in (14.461) . we obtain that 

||$i(t)|||<C 9 (l + log^) 6+Ar , 

||*i(t)ll!r < Cio2 i(f +1) (1 + log^) 2 (l + (1 + log*?)" max ||%_i(r)|| 



0<T<t 



Define now 



Zj-(f) = max «/||u(t)|| L2J , n) = max ||^-(t)||| 3 . 



- ^ y n«*v ;ilL»(n) - ^ II *i 

Then we have 



^i(t)<Cn(l + log£)( 6+A 0/ 4 , 
(*) < C^^+Da-'-- 1 ^ + log^^ 2 ) 2 ^" 1 max{l, ^(t)} . 
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In particular, putting yj(t) = max{l, Zj(t)} , we get 

(4.53) yi(t)<C n (l + log^)( 6+Ar )/ 4 , 

(4.54) y 3 {t) < (c 12 (l + lo gf? )(f +1 )) 2 J 2Hf+ 1 ) 2 - J % ._ 1 (t), for j > 2. 



Hence, passing to the logarithm in the inequality f)4.54p and summing up the result from 2 
to j , we obtain 

log%(t) < ^2-(log(c 12 (l + log^)^ +1 ) + l - (y + l) log 2) +bg(j/i(t)) 

i=2 ^ ' 

= log (c 12 (i + log *>) f +1 ) J2 2 ~ l + (f + lo s 2 E z2_l 

i=2 ^ ' i=2 

/ (6+iV) \ 

+ log(C r u (1 + logg) 4 J 
< log ( 2 (t + log ( Cm (1 + log ^ ) 

independently of j and t > 0. Hence, we deduce 

yj (t) < 2(t (c 14 (l + log,f +2 ), 

independently of j and t > 0. Choosing a proper (7, which is independent of g, we can 
conclude that 

sup J\\u(t)\\ LV < C(l + \ogg) 2+N . 

*>0,jGN V v 1 

Formula (I4.43P now immediately follows. 

4.4 Passage to the limit as n — > 00 

Our aim now is to pass to the limit in (I4.4|) — (I4.6p as n — > 00. 

From 1Q2]) . IjOijl . f l4~T5|) . IjOIj) . (IQ2| . it follows that, up to the extraction of some 
subsequence of n as n — >■ 00, there exist three functions u, 9 : (0, T) — > H , x : (0,T) — )■ H, 
such that we have (as a consequence of the generalized Ascoli theorem, see, e.g., [23J Cor. 8, 



p. 90]) 






(4.55) 


u n -> u 


weak star in ii^O, T; V 7 ) n L°°(0, T; if) n L 2 (0, T; V) 






and strongly in C°([0, T\\ V) D L 2 (0, T; ii), 


(4.56) 


6 n — >• 6 


weak star in L°°(0, T; if) n L 2 (0, T; V), 


(4.57) 


Xn^X 


weak star in L°°(0,T; V), 


(4.58) 




weak star in L°°(fi x (0, T)) d , 


(4.59) 




strongly in C°([0,T];H). 



29 



Moreover, uniformly bounded, it is easy to see that 

(4.60) Xn — > X weak star in L°°(0, T; V) and strongly in L 9 (Q T ), 

for every q G (1, oo), as n — > oo. Note that (cf. (14. 7p ) u > and 9 > a.e. Then it turns 
out that, at least for a subsequence, u n — >■ m a.e. Now, we denote by • , x) the inverse 
function of e with respect to the first variable, that is, e(i/}(w,x)) = w for all (w,x) £ 
[0, oo) x V(tp). Since e is continuous, increasing in 9, and such that e(6,x) > ejfl — e 2 for 
some constants e\ , e 2 > , we infer that ^ is continuous with a linear growth in [0, oo) xT>((p) 
and ip(0,x) = 0. The Nemytskii operator is therefore continuous in L 2 (Qt), and so this 
function is continuous and increasing with a linear growth. Hence, we have that 

9 n = tj) (u n - 9 n /n, Xn) -> ip(u, x) strongly in L 2 (0, T; H) . 

Hence, 9 = ip(u,x), or, equivalently, u = e(9,x)- Finally, we check that also 9 n converge 
strongly to 9 in L 2 (0,T; H) , at least for a subsequence of n — > oo. Indeed, from the 
definition of 9 n , we get 



(4.61) 



T 



o Jn 



E 

j=i "v-i 



9 n (x,t) - 9 n (x,t)\ 2 dxdt = 

r T/n " r i** 

/ \\9n(t)-9 , n f H dt+(^) £ / / 

JO XI 7 J=2 J(i Jtj-i 



9 n (x,t) - 9 n (x,t)\ 2 dtdx 



'3-1 



U-2 



,Xx,t) - 9 n (x,r)) dr 



dt dx . 



By (I4.56p . the first integral on the right-hand side of (I4.6ip can be estimated as follows: 

r-T/n q 
\ \\9n(t)- 9 , n \\ 2 H dt< -. 

Jo n 
Regarding the second term in (14.611) . we proceed as follows: 



J=2 
n 



tj-i 



U-2 



? n (x,t) - 9 n (x,r)) dr 



dt dx 



< 



n sr^ 



y-i 



\\9 n (t)-9 n (r)\\ 2 H drdt 



j=2 " Jt i-2 



< 



< 



n 



E 



t—U-a 



\9 n {t)-9 n {t-h)\\ 2 H dhdt 



j=2 Jt i-i J t—tj-i 
r2T/n / r T 



T Jo U ll^-^-^I^^J dh 



where we have used the new variable h = t — t . The last integral tends to , because 9 n 
converge strongly to 9, which is mean continuous in L 2 (Qt). This implies that 

9 n -> 9 strongly in L 2 (0, T; H), 
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and this allows us to pass to the limit in (j4.4|) — (I4.6p as n — > oo . Notice moreover that from 
estimate f l4.3(jp we also deduce that there exists a function ( : (0, T) — > V such that 

log# n -> C weakly in L 2 (0,T; V). 

Using the strong convergence of 9 n and the maximal monotonicity of the extended log graph, 
we also deduce that ( = log 9 G L 2 (0,T; V) , hence 9 > a.e. 



4.5 Conclusion of the existence proof 

Let us now introduce the limit problem obtained by passing to the limit as n — > oo in 
(I4.4p - fl4.6p in the previous subsection. 

PROBLEM (P)^. For fixed T > 0, and q > 0, find two functions 9 G L 2 (0,T;V) and 
X G L°°(tt x (0,T)) d , X t e L°°(tt x (0,T)) d such that, for t G (0,T), we have 

(4.62) ((e(9(t), X (t)))t,z)+ [ k(9(t), x(W9(t) ■ Vzdx + [ j(9(t) - 9 T ) zdA 

Jn Jan 

= - [ (\'(x)(t)xt(t) + /3Mx(t))) t -b[x}(t)xt(t)) zdx VzeV, 
Jn 

(4.63) Um)Xt{t) + 08 + 9(t))d<p( X )(t) 3 -X'( X )(t) - 0{t)a'{x){t) 
~ b[x\(t) ~ e x (9{t), X(t)) + 0(tX{0(t), X(t)) a.e. in Q , 

and for t = the functions e(9, x) and x satisfy the initial conditions 

(4.64) e(9,x)(0)=u , X (0) = X o ■ 



We thus have proved the following result. 

Proposition 4.2. Let Hypothesis \2.1\ hold true and let q > 0, fl g and s£ be defined as at 
the beginning of Section 14. J 1 Then Problem (P)^ has a solution in (0,T). 

In order to conclude the proof of Theorem \2.2\ we only need to remove the truncation 
q. This can be done using the Moser estimate f)4.43p . 

This bound for 9, indeed, allows us, for a suitable g > 1 satisfying 

C*(l + \og y +2N <q/2, 

with, e.g., C* = C 2 , to remove the truncation q and to conclude that the solution to 
Problem (P e ) is in fact a solution to Problem (P). 

Finally, let us note that we have now found a solution on (0, T) and we deduce that x 
is weakly continuous with values in V, e(9, x) £ C°([0, T];H), and 9 is bounded uniformly 
in time in L°°(n) due to the Moser estimate. Hence, we can continue the solution starting 
from time T and extend it on the whole time interval (0, oo) . 

This concludes the proof of Theorem 12.21 ■ 
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5 Uniqueness 

In this section, we prove Theorem 12. 51 Hence, we assume Hypothesis 12.31 First let us note 
that the lower bound for 9 (I2.18P directly follows by passing to the limit as n tends to oo in 
f)4.37p . Now we are ready to proceed with the proof of uniqueness and finally we will prove 
the regularity result (I2.2ip under the further assumption (12. 20 p . 



Uniqueness. In what follows, we denote by Ro, R\, R2, ■ ■ ■ suitable constants that possibly 
depend on T, but not on the solutions. We start with the proof of uniqueness of solutions. 
Equation (12. 2p is for (almost) all x G Q of the form (I2.10p . with 

a{ 9) - M 



(3 + 9 



g = £[9, X ] = -j^ (A'(x) + 9a' ( X ) + b[ X ] + e x (9, X ) ~ 9s x (9, X )) • 

Within the range < 9 < 9 and x £ ^ > c(v 9 ) > \Xt\ < C , of admissible values for the solutions, 
and, thanks to Hypothesis 12.11 (ii) (iii) (v) (vi), all nonlinearities in (I2.ip - (l2.2p are Lipschitz 
continuous. Using the notation from Theorem 12 .5[ we obtain, as a consequence of (12.171) . for 
a. e. (2, t) G Qoo the estimate 

(5.1) / \xt(x,r)\dr+\x(x,t)\ 

Jo 

< R (T) (\xo{x)\ + J (\9(x,t)\ + \x(x,r)\ + J \ X (y,r)\ dyj dr 
with some constant Rq(T) . Integrating over Q, and by Gronwall's argument, we obtain that 



(5.2) / \ X (y,t)\dy<R 1 / \ X o(y)\dy+ / / \6(y, r)\ dy dr 

Jn \Jn Jo Jn 

and hence, we get 



(5.3) / \xt&T)\dT+\x(x,t)\<R 2 (\xo(x)\+ / \9(x,r)\dr+ / \6(y, r)\ dy dr 
Jo \ Jo Jo Jn 

for a. e. x G Q and every t G [0, T] . In particular, 

(5.4) / / \Xt(x,r)\ dx dr < R 3 ( [ \x (x)\ dx + [ [ \9(x,r)\ dx dr) . 
Jo Jn \Jn Jo Jn J 

We now multiply (15. 3p by |x(x,t)| and integrate over Q to obtain for all t G [0,T] that 

(5.5) / \ X (x,t)\ 2 dx<Rj\\xo\\n+ f [ |^(x, r)| 2 dx dr) . 
Jn V Jo Jn J 

The crucial point is to exploit Eq. (12. ip properly. Notice first that we have 

(5.6) b[ X \xt(x, t) = 2B[ X ]t{x, t) + 2 / k(x, y) G'{ X (x, t) - X (y, t)) X t{y, t) dy . 
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We integrate the difference of the two equations fl2.ll) . written for (# l7 xi) an d (9 2 , X2) , from 
to t, rewriting the terms b[xi](Xi)t according to (15. 6p . Take z = K{6i) — K{9 2 ) in the 
resulting equation, where K{u) = j Q u k(s) ds, u £ M., and integrate it again over (0,t). 
Using the lower bound for 9, the Lipschitz continuity of all nonlinearities (<p is Lipschitz 
continuous on T>c((p) with constant C), the properties of K (cf. Hypothesis I2.1f iv)). and 
denoting 

Q(x,t) = [ (K(9 1 -K(9 2 ))(x,r)dr, 
Jo 

using (I4.37p . we obtain for each t £ (0,T) that 

t r 1 
2 - 



(5.7) k c v (w) / / \6(x,t)\ 2 dx dr + - / |V6(:r, t)| 2 dx 

<^(||0o||1t + IIxo||h + / / Ix0£,r)| 2 dxdr 
v Jo Jn 

+ [ [ [ I K{x,y)\xt{y,r)\\9{x,t)\dx dy drdC 
Jo Jo Jn Jn 

The last term on the right-hand side of the above inequality can be estimated, using (15 .4p . by 




*(x,y)\Xt(y,T)\\9(x,t)\ dx dy dr < R 6 / \9(x,t)\dx / / \xt(y, t)\ dy dr 
JnJn Jn Jo Jn 

<R 7 ( [ \9(x,t)\ 2 dx] ' ( [ \ X o(x)\ 2 dx + / / |#(x,r)| 2 dx dr) . 
\Jn J \Jn Jo Jn J 

Combining the last two inequalities again with the Gronwall's lemma, we obtain for each 
t £ [0,T] the estimate 

(5.8) / [ \9{x,r)\ 2 dx dr+ [ \Ve{x,t)\ 2 dx 



Jn 



< R% (\\6 \\h + Uo\\h + J q f n \^ X > T )\ 2 dxdT ) • 

We now multiply (15. 8p by 2_R 4 , add the result to (15. 5p . and see that Gronwall's argument 
can be applied again to arrive at the final estimate 

(5.9) / \x(x,t)\ 2 dx+ [ [ \9(x,r)\ 2 dxdr<R 8 (\\9 \\ 2 H + \\x \\n)- 

Jn Jo Jn v J 



Regularity. We prove now the regularity (I2.2ip for 9 under the further assumption (I2.20p . 
In order to do that, let us consider, instead of (I4.4p . the following approximated equation 

(5.10) J d t (\{t) + e{9 n {t), X n{t))j z dx + J k{9 n {t))V9 n {t) ■ Vz dx 

= ~ f ((A'(Xn)(t) +b\Xn}(t))d tX n(t) + 0d t (<f(Xn(t)))) Z dx . 
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Since (cf. Hypothesis I2.3( i)) the heat conductivity k is independent of x, we can now test 
f)5.10p by K(9 n )t : where K(6 n ) = J Q n k(s) ds and, integrating over (0,t), we obtain, using 
f)4.37p . the estimate 

ko + c v (w^j J \{d n ) t \ 2 dxdr + \\\V{K{9 n )){t)f H 

< J||V(A(M)lll + C^i / / \(e n )t\dxdr. 
2 Jo Jn 

Here C is a bound in L°°(Q T ) for the term (A'(xn) + b[x n ])dtXn + /3d t (^(Xn)) that we 
already obtained in Theorem I2.2| and ko , k\ are the positive constants introduced in Hy- 
pothesis ET^iv). Using now assumption (12. 20 j) and the conclusions of Theorem l2.2| we obtain 
the (independent of n) bound 

||(#n)t||z, 2 (0,T;H) + || K{6 n ) \\ L oo ( ,T;V) < C , 

which leads immediately (due to Hypothesis I2.1( iv)) to the desired estimate (12.21 j) . 

With this, Theorem 12.51 is proved. ■ 

Remark 5.1. It seems to the authors that, in the case when k depends only on 6, if we 
replace the boundary condition (II. 3p with 

(5.11) k(6) W • n + j(K(6) - K(6 V )) = onffi, 

where K is the primitive of k as in this section, then the uniqueness result holds for a 
general nonnegative 7 . Moreover, the regularity estimates can be adapted in order to cover 
also this case: the details of this proof will be developed in a subsequent paper. However, 
let us point out that the regularity properties ( 12 . 2 1 [) remain true for more general 7 also in 
the case of boundary condition (II. 3p : the argument relies on a suitable use of the Gronwall 
lemma. 
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